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Abstract. We study electroweakly interacting particles in rotating matter. The existence 
of the electric current along the axis of the matter rotation is predicted in this system. This 
new galvano-rotational effect is caused by the parity violating interaction between massless 
charged particles in the rotating matter. We start with the exact solution of the Dirac equa¬ 
tion for a fermion involved in the electroweak interaction in the rotating frame. This equation 
includes the noninertial effects. Then, using the obtained solution, we derive the induced 
electric current which turns out to flow along the rotation axis. We study the possibility 
of the appearance of the galvano-rotational effect in dense matter of compact astrophysical 
objects. The particular example of neutron and hypothetical quark stars is discussed. It 
is shown that, using this effect, one can expect the generation of toroidal magnetic fields 
comparable with poloidal ones in old millisecond pulsars. We also briefly discuss the gen¬ 
eration of the magnetic helicity in these stars. Finally we analyze the possibility to apply 
the galvano-rotational effect for the description of the asymmetric neutrino emission from a 
neutron star to explain pulsars kicks. 
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1 Introduction 

The importance of noninertial effects for various areas in modern physics cannot be under¬ 
estimated. Some of the examples of these effects are mentioned in ref. [1]. One of the most 
common manifestations of noninertial effects is the description of physical processes in a ro¬ 
tating frame. One can expect the appearance of additional interesting phenomena if, besides 
the matter rotation, there is a parity violating interaction in the system. In the present work 
we will show that, in this situation, an electric current flowing along the rotation axis can be 
induced. We shall call this phenomenon as the new galvano-rotational effect (GRE). 

Previously the generation of an electric current due to nontrivial topological effects was 
studied mainly in connection to the chiral magnetic effect (CME) [2]. In that case a nonzero 
current can be induced in the system of massless particles embedded in an external magnetic 
field [3], provided there is an imbalance between left and right particles. Recently, in ref. [4], 
we showed that the electric current can be generated even at zero chiral imbalance if charged 
particles are involved in the parity violating interaction. 

We will apply the new GRE in astrophysical media to generate a toroidal magnetic 
field (TMF) inside a compact star. Although stellar TMFs cannot be observed directly, they 
are an internal ingredient of various astrophysical objects. For example, the most plausible 
explanation of the 22 yr solar cycle is the oscillation between poloidal and toroidal components 
of the solar magnetic field [5]. Moreover, purely poloidal or toroidal stellar magnetic fields 
were shown in refs. [6, 7] to be unstable. Thus TMF is inherent to a magnetized star. In this 
work we show how TMF can be generated in an old millisecond pulsar basing on GRE. 

Since macroscopic fluxes of electroweakly interacting particles are produced by GRE, 
we can try to apply this effect to explain linear velocities of millisecond pulsars. It is known 
from astronomical observations [8] that pulsars possess great linear velocities. Nevertheless 
physical processes underlaying pulsar kicks are still unclear. It might be reasonable to use 
GRE to account for pulsars kicks due to, e.g., anisotropic neutrino emission since linear 
velocities of pulsars were reported in ref. [9] to be correlated with their angular velocities. 
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This paper is organized in the following way. First, in section 2, we briefly describe 
the Standard Model interaction between leptons and quarks in flat space-time. Then, in 
section 3, we derive the Dirac equation for a fermion which interacts electroweakly with a 
rotating background matter. For this purpose we write down this Dirac equation in the 
corotating frame, using the method of an effective curved space-time. The exact solution of 
the Dirac equation for an ultrarelativistic fermion, accounting for the noninertial effects, is 
obtained in section 3. In section 4, we establish GRE, which, in this situation, consists in 
the appearance of the electric current along the rotation axis. We calculate this current in 
section 4 using the exact solution of the Dirac equation obtained in section 3. In section 5, we 
apply GRE for the generation of TMF and the magnetic helicity in compact rapidly rotating 
stars. Finally, in section 6 , we try to use GRE to produce anisotropic neutrino emission from 
pulsars to explain their great linear velocities. 

In section 7, we summarize our results and compare them with the findings of other 
authors. Some details of the derivation of the electric current in rotating matter are provided 
in appendix A. 


2 Electroweak interaction of fermions in flat space-time 


In this section we shall briefly remind the description of the electroweak interaction between 
leptons and quarks in flat space-time. 

Let us consider a medium consisting of electrons, electron neutrinos v e as well as u and 
d quarks. We shall assume that quarks are both in confined states, forming nucleons, and 
hypothetical free particles. The effective Lagrangian for the electroweak interaction in this 
system in the Fermi approximation has form [10], 


-Ccff 


4G f 

75 


j/y/t + K^Rp, , 


( 2 . 1 ) 


where Gp ~ 1.17 x 10 -5 GeV -2 is the Fermi constant, J^ is the charged current, and AT is 
the neutral current. 

In the considered system of elementary particles, the charged current has the form, 


Jfl — Vud'4 ) u'YlJ,'4 , d + VVeT/i Vt 


( 2 . 2 ) 


where ip u .d are the wave functions of u and d quarks, tpv e ,e are the wave functions of v e and 
an electron, 7 ^ = 7 ^ (l — 7 s ) /2, 7 ^ = ( 7 °, 7 ) are the Dirac matrices, 7 s = i 7 ° 7 1 7 2 7 3 , and 
Vud ~ 0.97 is the element of the Cabbibo-Kobayashi-Maskawa matrix. The neutral current 
can be expressed as 

K v = ^ [ € /V’/7m^/ + e/^/7 *Vt] > ( 2 - 3 ) 
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where 7 ^ = 7 ^ (l + 7 5 ) /2 and 
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Here £ = sin 2 ~ 0.23 and 0w is the Weinberg angle. 


(2.4) 
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Table 1 . The values of the effective potentials Vl.r in eq. (2.5) for various channels of the scattering 
of a test fermion off background particles. Here n e is the electron density, n njP are the densities of 
neutrons and protons, and n Ut d are the densities of u and d quarks. 


In this section we shall consider the case when background fermions are at rest and 
unpolarized. While averaging the currents in eqs. (2.2) and (2.3) over the Fermi-Dirac dis¬ 
tributions (•••), in this case we get that only = n f 7^ 0, where rif is the invariant 

number density of these fermions. The quantities (V’/T^Z’/) = 0 and ('*A/7 /i 7 5 V ; /) = 0 are 
vanishing since they are proportional to the macroscopic velocity and the polarization. 

After averaging over the ensemble of background particles, we can rewrite eq. (2.1) in 
the form, 

Aff = -i> [tS'Vl + 7o R Vr] ^ (2.5) 

where ^ is the wave function of a test fermion which undergoes a scattering off background 
particles and the effective potentials Vl,r are given in table 1 for any scattering channels. 

When we consider the electron and v e scattering off nucleons, u and d quarks are 
confined inside neutrons or protons. In the case of electron-nucleons interaction, only the 
neutral current contributes to eq. (2.5). If we study the scattering u quarks off d quarks, and 
vice versa, as well as the v e scattering off background fermions, both the charged and the 
neutral currents give the contributions to eq. (2.5). To derive the expression for Vl for (ud), 
( du ), and (z^e) interactions on the basis of eq. (2.2), we use the Fierz transformation. We 
also note that we consider the v e interaction with electroneutral matter where n e = n p . 

3 Dirac equation for a fermion interacting with rotating matter 

In this section we shall find the exact solution of the Dirac equation for a fermion interacting 
with a rotating matter by means of the electroweak forces. The obtained solution includes 
noninertial effects. 

Let us discuss the interaction of a fermion with matter rotating with the constant 
angular velocity oj. In this case we cannot directly apply the results of section 2 taking 
{^}fyi /)/) ~ vj / 0, where vj = (u> x r) is the fermions velocity, while averaging over the 
ensemble of background particles. In the situation, when matter moves with an acceleration, 
one should account for possible noninertial effects. 

Nevertheless we can still choose a noninertial reference frame where matter is at rest. 
Assuming that matter is unpolarized, we get that only ^ 0 in this reference frame. 

Thus, formally we can use the effective potentials derived in section 2. For the first time this 
approach was put forward in ref. [11], where the neutrino interaction with a rotating matter 
was studied. 
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It is known that the description of a particle in an accelerated frame is analogous 
to the motion of this particle in the curved space-time or the interaction with an effective 
gravitational held. For example, when we study the motion in the rotating frame, the interval 
takes the form [ 12 ], 

ds 2 = g^dx^dx 1 ' = (1 — co 2 r 2 )dt 2 — dr 2 — 2ur 2 dtd(p — r 2 d(f> 2 — dz 2 , (3.1) 

where g^ v is the metric tensor of the effective gravitational held. Here we use the cylindrical 
coordinates x^ = ( t,r,4>,z ). 

Using eq. (2.5), we get that the Dirac equation for a test fermion, with the mass m, 
involved in the parity violating interaction and moving in a curved space-time, has the form 
(see also ref. [ 11 ]), 

f v , J ' V M I 

[i7 /i (x)V /i - m] ij> = 7^(x) | -y [l ~ 7 5 (^)] + -y [l + 7 5 (u>] j (3-2) 

where 7 /i (x) are the coordinate dependent Dirac matrices, = d^ + T^ is the covariant 
derivative, T M is the spin connection, 7 5 (x) = — (i/dl)E fJlUa P'y fl {x)'y v (x)'y ot {x)'yp(x), E^ a P = 
£ [ivctpj_g ^j ie C ovariant antisymmetric tensor in curved space-time, g = det(g fiu ), and 
U/ R are the effective potentials. Note that, since we choose a corotating frame, then D r ° R = 
Vl,r / 0 and R = 0, where Vl,r are given in table 1. Analogous Dirac equation was also 
discussed in ref. [13]. 

One can check that, using the following vierbein vectors: 


e 0 M = (1,0, —uj, 0), e/ = (0,l,0,0), e/ = (0,0,l/r,0), e 3 M = (0,0,0,1), (3.3) 

the metric tensor in eq. (3.1) can be diagonalized, g a b = e^e^g^, where r] a f) = diag(l, —1, 
— 1,-1) is the metric in a locally Minkowskian frame. 

Let us introduce the constant Dirac matrices in a locally Minkowskian frame by j a = 
e °ijH IJ (x)- As shown in ref. [11], 7 5 (x) = ry ^ 1 ^^ 3 = 7 s does not depend on coordinates. 
The spin connection in the Dirac eq. (3.2) has the form [14], 


r " = -4 


® ^a6/x? ^ab/i 


(3.4) 


where a a b = (i/2) [ 7 ^, 75 ] — are the generators of the Lorentz transformations in a locally 
Minkowskian frame and the semicolon stays for the covariant derivative. The explicit calcu¬ 
lation on the basis of eq. (3.4) shows that the nonzero components of the connection one-form 
w s (, = %dak are 

wi 2 M = -w 2 i m = (w,0,1,0). (3.5) 

Using eqs. (3.4) and (3.5) we get that ry^/a/T^ = /2r. 

Using the definition of 7 “, the Dirac eq. (3.2) takes the form, 


= 

V = 


( 7 ° - wry 2 ) (vy - Ua 7 S ) V’, 

i 7 ° (do - ud^) + ky 1 + y j 


+ 17^ + 17% 
r 



(3.6) 


where Uv,A = (Ul ± Ur) /2 are the vector and axial parts of the effective potentials. Note 
that the operator V in eq. (3.6) is analogous to that recently studied in ref. [15]. 
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Since eq. (3.6) does not explicitly depend on t, (j), and z, we shall look for its solution 
in the form, 

if} = exp (—iEt + iJ z (j) ± ip z z) i/j r , (3-7) 

where ip r = i^ r { r ) is the wave function depending on the radial coordinate. The values of J z 
in eq. (3.7) were found in ref. [16] to be ±1/2, ±3/2,.... 

It is convenient to rewrite eq. (3.6) as 


where Q a = q a 


[7 a Qa -m + V]lp r = 0, 

(feffTlgg, q e s is the effective electric charge, q a = (E ± J z uj 


^eff= ( 0 ,-^,— 
Qeff 


Vycor — — 
r 


(3.8) 

Vy,-id r , 0 ,p z ), 

(3.9) 


is the potential of the effective electromagnetic field, and V = Va ^y° — wry 2 j y 5 . 

Let us look for the solution of eq. (3.8) in the form, ip r = [7 a Q a ± m — 17] <h, where <!> 
is the new spinor. The equation for <h reads 


1 \ 


2 r J 


J z 


d r ± — I ± (E ± J Z UJ — hy) — I- Vyur ) ± I Vyw H—7 ) S3 


J z 


± 2Va7 5 


— — Vyur ) ur — (E ± J z u — Vy) ± ^S 2 - p" 
r J 2 

+ Va (l — cu 2 r 2 ) + 2mV\ ( 7 ° — wry 2 ') 7 ® — m 2 


<h = 0, (3.10) 


where S 3 = 7 ° 7 3 7 5 . 

The solution of eq. (3.10) can be found for ultrarelativistic particles. In the limit m —» 0, 
we can represent <h = vp in eq. (3.10), where p = p{r) is a scalar function and v is a constant 
spinor satisfying = crv and y 5 u = x v i with a = ±1 and x = ±1) since both S 3 and y 5 
now commute with the operator of eq. (3.10). 

Let us first study left particles, (1 ± y 5 )^ = 0, corresponding to x = ±1. The case of 
right particles with x = — 1 can be studied analogously. Using the new variable p = |VL|wr 2 
we can write the equation for p a as 


1 


4 P 


a 


1 


P 1 


pd + d p - — (l ± -sgn(U L ) - - ---- I - -sgn(U L ) - - ± k 


4 2 


a 


1 


Pa = 0, (3.11) 


where 


k = 


4|U l |w L 


E A ± 2 E ( J z u - V L ) -p'i-m' z + (U L ) Z ± JW ± V L ua 


a 


+ o [ l - oSgnCVL) - o 


(3.12) 


In eq. (3.11) we choose J z = (1/2 — ^sgn/Ii,), where l = 0, ±1, ±2,.... 

Assuming that p a —>• 0 at r —>• 00 , we get that p + = 7jv, s and p- = T/v-gs, for Vl > 0, 
as well as p + = In-i,s and p_ = In,s, for Vl < 0 . Here AT = 0 , 1 , 2 ... , s = N — l, and 
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In,s = In,s{p ) is the Laguerre function 1 . The energy spectrum can be found if we take that 
ac = N in eq. (3.12). We can present E is the form, 


E = V L + (l - ^ wsgn(U L ) ±S, £ = ^p 2 z + m 2 + 4\V L \Nuj. (3.13) 

One can see that the neutrino energy depends on the sign of Vl- Note that one should 
understand m ^ 0 in eq. (3.13) in the perturbative sense. 

The total radial wave function ip r can be found in the explicit form if we choose the 
spinors v a as 

^ = (1,0,0,0), vZ = (0,l,0,0). (3.14) 

In eq. (3.14) we assume that the Dirac matrices are in the chiral representation [17], 


7 


o 




7 


f 0 a k \ 

V -<Tk 0 ) ’ 


7 



(3.15) 


where are the Pauli matrices. 

Using eq. (3.14), we get the radial wave functions corresponding to different spin pro- 


jections, 

7- 

’’ H- 

II 

c: 

where 

n =( 

f Pz ~ 

and 




7+ = n 


<P+ 

0 


p z -E- J z uj + V L 


rj- = n 


— iR± 


0 

<P- 


-i R 


,= F 


Pz ~ E - J z u> + Vl 


R + =d r + ^ + y “ v ^ r = V\ v l\ up ^2 d p - - 1^ 


1 J. 


R- —d r + — —— + Vli vr — sj | Vi,|cap ( 2d p + — + 1 ) . 


In eq. (3.17) the upper signs stay for Vl > 0 and the lower ones for Dl < 0. 

The properly normalized two component spinors in eq. (3.16) have the form, 


(3.16) 

(3.17) 


(3.18) 


_ (-1 ( T [ £ T Pz\ I n.s \ ri = r ( 2i VW^In,s A 

\ T [£±Pz] In- 1,s) 

for Ul > 0 , and 

_n ( ~F Pz] In -i, s ^ f — 2i\/|ViJu;lVIjv_i jS \ 

V+ ° + V 2iy/\V[kjNI NiS J ’ V T [£± Pz \In,s ) 


(3.19) 


(3.20) 


for Vl < 0. The signs in eqs. (3.19) and (3.20) are correlated with the sign in eq. (3.13). 

It should be noted that the spinors 77 + and ??_ in eqs. (3.19) and (3.20) are linearly 
dependent as it should be for ultrarelativistic particles. Thus, one can use any independent 
pair of rj + and ??_. As usual, we shall attribute ??_ with the upper sign to a particle degree 
of freedom and r) + with the lower sign to antiparticles. This choice of independent spinors 

1 The definition of the Laguerre function is given, e.g., in ref. [11]. 
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Values of Vl,r 

Allowed value of p z 

Left particles 

V L >0 

Pz <0 

v L <0 

Pz> 0 

Left antiparticles 

V L >0 

Pz> 0 

v L <0 

Pz < 0 

Right particles 

Ur > 0 

Pz> 0 

Ur < 0 

Pz <0 

Right antiparticles 

Ur > 0 

Pz <0 

Ur < 0 

Pz> 0 


Table 2. Allowed values of p z in the ground state N = 0 for different values of the effective potentials 
Vl,r for left and right particles and antiparticles. 


is convenient for N > 0. If N = 0, one can better use 77 + for Vl > 0 and for Vl < 0 as 
independent degrees of freedom. 

Using the normalization condition for the total wave function, 

J ^N 1 s,pS x )^ N '’ s '’P'S x )'f-9 d3 x = $NN'5 ss '5(p z - p' z ), (3.21) 


which includes the dependence on p z and (j), we get the coefficients C a in eqs. (3.19) and (3.20) 
as 


,2 _ |Ui> 

CT 2tt£(£t crpz)' 


(3.22) 


which is valid for any sign of Vl. 

On the basis of eqs. (3.19) and (3.20) one can notice that, at N = 0, p z is correlated 
with the particle helicity. Using eq. (3.13) with m = 0 as well as eqs. (3.19) and (3.20), one 
can find the possible values of p z at IV = 0. For the convenience, they are listed in table 2. 
Note that at N > 0, —00 <p z < + 00 . 

It should be noted that p z in table 2 is a formal quantum number. The physical value 
of p z for antiparticles is opposite to that shown in table 2 : pi pliyb '' = —p z . Otherwise the 
electric charge would not be conserved. 

Right particles can be treated in the same way as left ones. That is why we just present 
only the final results. The expression for the energy has the same structure as eq. (3.13) with 
the replacement Vl —>• Ur- In eq. (3.16) one has 'tpf 1 = (£±, 0 ) T , where 


<Trv 

+ 

II 

-p 

f T [£ ± Pz] In,s \ 
v 2V|Ur|u;JVUv_i >s J 

p _ ( — 21i/|Vr \u>NI n ,s 

- V =F [S =F p z \ In-i,S 

for Ur > 0, and 



£+ = <?+ ( 

{ T [£ ±Pz\ In-i, s \ 
v -2V|Ur| ujNI Nt J’ 

p _ n ( 2iyy\Vn,\u}NlN-i,s 
~ \ T [£ T Pz] In, s 


(3.23) 


(3.24) 
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for Vr < 0. The argument of the Laguerre functions is p = |UR,|wr 2 now. The new normal¬ 
ization constant in eq. (3.23) and (3.24) reads 


,2 = IVrI^ 

a 2n£(£ ± crpz) 


(3.25) 


The signs in eqs. (3.23)-(3.25) are correlated with the signs in the expression for the energy. 

As in the case of left particles, for right fermions we have that — oo < p z < +oo at 
N > 0. The allowed values of p z at N = 0 are shown in table 2. 


4 Induced electric current along the rotation axis 

In this section we show that there is a nonzero induced electric current flowing along the 
rotation axis in the system of electroweakly interacting particles. In our calculation we shall 
use the exact solution of the Dirac equation obtained in section 3. 

In section 3 we already mentioned that there is a correlation between p z and the he- 
licity at N = 0. Thus one expects that there can be macroscopic fluxes of particles in the 
rotating matter. Let us first examine this issue for left fermions. We shall calculate the mean 
hydrodynamic currents of particles and antiparticles with respect to the coordinates x 11 in 
the rotating frame. These currents have the form, 

OO ,, 

3LfJ = / dp g ‘^f{x)i/;f^j(£±Vj J ), (4.1) 

N , s = 0 ' 

where p^j j{E) = {exp(/3[£'=p pi]) + 1} _1 is the Fermi-Dirac distribution for fermions, with 
the lower sign staying for antifermions, /3 = 1/T is the reciprocal temperature of the fermion 
gas, and pi, is the chemical potential of left particles. The spinor in eq. (4.1) corresponds 
to the exact solution of the Dirac equation in eqs. (3.19) and (3.20). Note that, for the first 
time, this method for the calculation of the current was proposed in ref. [3] . 

We will be interested in the expression for j^'j linear in uj. That is why we use £ + 
Vl instead of the total particle energy, cf. eq. (3.13), in the distribution function. The 
contribution of the noninertial part of the energy ~ o;sgn(l^ J] R)(Z — 1/2), see eq. (3.13), to 
the currents is computed in appendix A. We should study only j^j- since it is this component 
of the current that is linear in u. 

Using the orthogonality of Laguerre functions, 

OO 

^2In,s{p)In',s(p) = 5nn', (4.2) 

s =0 

as well as eqs. (3.3), (3.7), (3.15), (3.16), (3.19), (3.20), and (3.22) we get that 

w = ±^lhL, (4.3) 

n £ 7T 

s=0 

where the upper sign stays for particles and the lower one for antiparticles. 

On the basis of eq. (4.3) and table 2 we find that only the lowest level IV = 0 contributes 
to the current. Finally we get that 

j'tf = / dpptfip + V-L). (4.4) 

71 Jo 
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Analogously to eq. (4.4) we can obtain the following expression: 


J'l/ = “ ~~ J o d PP L f(P ~ VL), (4.5) 

which is valid for antifermions. 

The contribution to the hydrodynamic current from right femions is analogous to 


eqs. (4.4) and (4.5). It is 

Vo. 

£dco 

II 

Erw 

7T 

POO 

/ d PPf(P + VR), 

Jo 

(4.6) 

for particles, and 

J R / “ 

Fru; 

7r 

POO 

/ d ppf(p-V R ), 

Jo 

(4.7) 


for antiparticles. Here p^' j{E) are the distributions of right particles and antiparticles which 
can be obtained from j{E) by replacing /zl —> /zr, where pa is the chemical potential of 
right fermions. 

Now, using eqs. (4.4)-(4.7), we can obtain the expression for the third component of 
the electric current as J 3 = qf (j R f — j 3 j + j^f — j R jJ, where qf is the electric charge of 
the fermion / including the sign. For example, q e = —e for an electron, q u = 2e/3 for an 
u quark, and qd = —e/3 for a d quark. Here e > 0 is the absolute value of the elementary 
electric charge. In the expression for J 3 , we use the convention that the direction of the 
electric current coincides with the motion of the positive electric charge. Finally we get for 
the electric current, 

j _ ( VrI _ 1r _ yL/i L ), (4.8) 

7r 

where we restore vector notations. We remind that in eq. (4.8) we keep only the terms linear 
in id and Vl,r- Some details of the derivation of eq. (4.8) from eqs. (4.4)-(4.7) are provided 
in appendix A. 

We can attribute the existence of the induced electric current in rotating matter, where 
the parity violating interaction is present, to the new GRE; cf. section 1. 

5 Generation of TMF and magnetic helicity in a pulsar 

In this section we apply GRE for the calculation of TMF inside a pulsar. We also briefly 
consider the generation of the magnetic helicity in a compact rotating star. 

If we consider a rapidly rotating compact astrophysical object, like a neutron star (NS) 
or even a hypothetical quark star (QS), then the mechanism described in section 4 will induce 
the electric current along the rotation axis of such a star. We shall suppose that this current 
forms a closed circuit connected somewhere at the stellar surface. Thus, using the Maxwell 
equation (V x B) = J, we get that this current should induce a TMF, Bt OI ~ RJ, where 
R ~ 10 km is the stellar radius. 

It should be noted that a compact star typically has a poloidal magnetic field B po i, 
which is measured in astronomical observations. For instance, the radiation of a pulsar 
can be explained by the emission of electromagnetic waves by the rotating magnetic dipole 
associated with R po i, provided there is a nonzero angle between B po i and id. However, as 
shown in ref. [18], using general arguments for the magnetohydrodynamic equilibrium of an 
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axisymmetric NS, a purely poloidal magnetic field configuration turns out to be unstable. 
Thus an internal nonzero TMF should exist in a compact star. 

Let us first consider the generation of TMF in NS composed of degenerate electrons 
and nucleons, like neutrons and protons. In this case u and d quarks are confined inside 
nucleons. The typical electron density in NS is n e ~ 9 x 10 36 cm' 3 , which corresponds to the 
electron fraction Y e ~ 0.05. It gives the chemical potential of electrons fi e ~ 125 MeV m e . 
Therefore electrons are ultrarelativistic whereas neutrons and protons are nonrelativistic. 
Note that the nonzero electron mass was shown in ref. [19] to slightly contribute to the 
induced electric current. Thus we can assume that electrons are approximately massless in 
NS and the results of sections 3 and 4 are valid. 

Since the chiral symmetry is unbroken, we can consider left and right chiral projections 
as independent degrees of freedom. For simplicity we shall take that left and right electrons 
are in equilibrium with m, ~ //r ~ fi e . The situations, when the chiral imbalance = 
(fiR — /^l)/2 / 0 is important in NS, are studied, e.g., in refs. [4, 19]. Since n v <C n n in NS, 
(Vi, — Ur) = Gprin/y/2 ~ 12 eV for n n « 1.8 x 10 38 cm" 3 . Eventually, taking that u ~ 10 3 s" 1 
and R ~ 10 km as well as using eq. (4.8), we get that B tor ~ 2.5 x 10 8 G can be generated in 
a rotating NS. 

The obtained value of Bt or is comparable with B po i < (10 8 —10 9 ) G in weakly magnetized 
old millisecond pulsars [20]. It should be noted that the stability of magnetic fields in NS 
can be reached if 0 < E tor / E mag < 0.07 [18], where E toI ~ L?/ or is the energy of TMF 
and E mag (B^ or + -Bp ol ) is the total magnetic energy. Our estimate for B tor satisfies this 
criterion. 

Let us discuss the creation of TMF in a hypothetical QS. Although QSs have not been 
observed yet, their properties are actively studied theoretically [21]. Various models of QS 
predict that it consists of free u and d quarks with some admixture of s quarks. After the 
analysis of various equations of state of QS matter, the strangeness fraction was found in 
ref. [22] not to exceed ~ 0.3. Thus we can approximately omit the contribution of s quarks 
in the calculation of the induced current. 

We shall suppose that inside QS we have n u = no/3 ~ 0.6 x 10 38 cm -3 and n^ = 2no/3 ~ 
1.2 x 10 38 cm~ 3 , where no ~ 1.8 x 10 38 cm -3 is the nuclear density. At such high densities 
the chiral symmetry can be unbroken in QS [23], allowing one to consider the independent 
chiral projections of the wave functions of u and d quarks. Therefore we can again use the 
results of sections 3 and 4. As in case of NS, we can also assume that left and right quarks 
are in equilibrium, fi^ d ~ fi^ d ~ li u ,ch just for simplicity. Here fi u = (1/3) 1 / 3 /Uo ~ 235 MeV 
and fid = (2/3) 1//3 /io ~ 296 MeV, where fio = 339 MeV. 

Finally, using, in eq. (4.8), the adopted values of densities and chemical potentials of 
quarks, the values of Ll,r in table 1 as well as for ui 10 3 s 1 and R. 10 km, we get that 
-Btor ~ 4.9 x 10 8 G can exist inside a rotating QS. The obtained value of TMF is slightly 
greater than L> t0 r for NS. Note that the derived strength of TMF is also in agreement with 
stability criterion obtained in ref. [18] for old weakly magnetized millisecond pulsars [20]. 

It should be noted that, in our estimate of TMF in QS, we account for only (ud) and 
(du) contributions to the electric current. Using the analogy of the rotating electroweak 
matter with the presence of an effective magnetic field [24] (see also eqs. (3.8) and (3.9) in 
section 3) and the results of ref. [25] , we get that (uu) and (dd) interactions do not contribute 
to the current in eq. (4.8). 

We should mention that, in the generation of TMF in a compact star, we discuss a 
thermally relaxed stage in the evolution of this astrophysical object. This approximation is 
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valid since we consider old millisecond pulsars with ages (10 8 — 10 10 ) yr [20]. It means that 
we discard any possible effects related to turbulence which should be treated on the basis of 
the Navier-Stokes equation. In our analysis we also do not consider a differential rotation 
either. 

The creation of TMF in a compact star is closely related to the problem of the generation 
of the magnetic helicity defined as 

H = J d 3 z(A-B), (5.1) 

where A is the 3D vector potential. If there is configuration of magnetic fields in a star con¬ 
sisting of toroidal and poloidal fields, then H in eq. (5.1) has the form [26], H = 2L<f> t or < I ) poi, 
where <h t0 r aR d ‘hpol are the fluxes of toroidal and poloidal fields and |L| = 1 is the linkage 
number. The magnetic helicity is a conserved quantity in a perfectly conducting medium. 
It is this fact which provides the stability of a poloidal field in a compact star. Note that 
another mechanism for the generation of the magnetic helicity in a nonrotating NS, based 
on the electron-nucleon electroweak interaction, was recently proposed in refs. [4, 27]. 

6 Pulsar kicks due to the asymmetric neutrino emission 

In this section we shall use GRE for the description of the asymmetry in the neutrino emission 
from NS. We shall also consider the applicability of our results to explain linear velocities of 
pulsars. 

It is well established that some pulsars have great linear velocities up to 10 3 km-s 1 [8]. 
There are various models for the explanation of this phenomenon based on, e.g., the asym¬ 
metric electromagnetic radiation [28] and the asymmetric explosion leading to the anisotropic 
neutrino emission [29]. We also mention ref. [30], where the asymmetry in neutrino oscil¬ 
lations in matter and an external magnetic field was used to account for pulsar kicks. The 
idea that anisotropically emitted electrons, owing to CME, pass the momentum to NS, was 
discussed in ref. [19]. Nevertheless the origin of peculiar velocities of pulsars is still unclear. 

In ref. [9] it was established that linear velocities of pulsars are correlated with their 
angular velocities. Thus we may try to apply GRE, which predicts particle fluxes along 
the rotation axis, to explain pulsar kicks. However, unlike ref. [19], we shall examine the 
possibility of asymmetric neutrino emission since it is not very clear how charged particles 
can escape NS. 

Using eqs. (4.4) and (4.5), we get the total hydrodynamic current of neutrinos along 
the rotation axis as 

ujT r°° ( 1 i \ 

jl = GFHn W2 F{y)l F{y) = .lo (^Ti + ^TT) dx = 2ln{1 + ey) - y ’ (fU) 

where y = (/tl ~ Vl)/T and T is the neutrino temperature. To derive eq. (6.1) we take into 
account that Ul ~ —Gprin/y /2 for v e (see table 1) and n n n e . Note that both /tr and Vr, 
are equal to zero for ultrarelativistic neutrinos. That is why we account for the contribution 
of only left neutrinos in eq. (6.1). 

The more complicated structure of the hydrodynamic current in eq. (6.1) compared to 
that of the electric current in eq. (4.8) is owing to the fact that j 3 ^ has the same direction 
as j'lj, whereas J 3 j is directed oppositely to Jf j-. The function F(y ) is plotted in figure 1. 
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Figure 1. The function F in eq. (6.1) versus y. 

It is interesting to mention that -F(O) ~ 1.4. Thus, there is a nonzero neutrino flux even at 
/iL = 0. 

The simulations carried out in ref. [31] show that there is a significant nonzero neutrino 
asymmetry A n Ve = n Ve — np e , owing to direct Urea processes, which lasts up to t\ ~ 0.1s 
after the onset of the supernova collapse. When 0 < t < t\, the typical neutrino energy 
E v ~ lOMeV and chemical potential of v e is /iL ~ lOMeV [32]. At t\ < t < t 2 ~ 10 2 yr other 
neutrino species start to be emitted resulting in the diminishing of /iL- When t 2 < t < £3 ~ 
10 6 yr only vv pairs can be emitted in modified Urea processes [33] leading to //l = 0. For 
simplicity we shall assume that T 10 8 K = const and E u ~ lOkeV [33] at this stage of the 
NS evolution. 

The total momentum carried away by neutrinos during the time interval At is P ~ 
j 3 E u SAt, where S = nR 2 is area of the equatorial cross section of NS and and R ~ 10 km 
is the NS radius. Using eq. (6.1), one gets that the greatest P is achieved at t 2 < t < f 3 - 
As a result, NS will get a recoil velocity v = P/M , where M is the NS mass. We shall take 
M = 1.44Mq, where Mq « 2x 10 33 g is the solar mass, to get the upper bound for v. Taking 
that w 10 3 s _1 , we get that v ~ 10 16 cm • s 1 . The obtained value is sure to be beyond 
the possibility of modern astronomical observations. Therefore, despite there is a nonzero 
anisotropy in the neutrino emission in NS owing to GRE, this effect is unlikely to result in 
any testable phenomena. 

7 Discussion 

In conclusion we note that in the present work we have studied the evolution of particles, 
involved in the parity violating electroweak interaction, in the rotating matter. In section 3, 
we have obtained the new exact solution of the Dirac equation for a test ultrarelativistic 
particle, which account for the noninertial effects. Then, in section 4, we have computed 
the induced electric current along the rotation axis on the basis of the exact solution of the 
Dirac equation. In section 5, we have applied our results for the generation of TMF and 
the magnetic helicity in compact rotating stars. Finally, in section 6 , we have considered the 
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production of the anisotropy in the neutrino emission from NS, owing to GRE, and examined 
the applicability of this effect for the explanation of peculiar velocities of pulsars. 

Several new results have been obtained in this work. Firstly, we mention that the 
vierbein vectors in eq. (3.3) have never been used previously in the Dirac eq. (3.2), which 
accounts for the electroweak interaction with background matter in curved space-time. An¬ 
other veirbein was recently used in ref. [11]. However, the choice of the vierbein in the present 
work is likely to be more appropriate for ultrarelativistic particles in a rotating frame. In 
particular, here we have obtained the correct form of the “centrifugal” energy, or the energy 
of the rotation—angular momentum coupling, E c { = — (J • u>): cf. eq. (3.13). The obtained 
expression for E c f coincides with the result of ref. [34] derived on the basis of the general 
analysis. The form of E c f obtained in ref. [11], where another vierbein was used, is slightly 
different. Therefore, the vierbein adopted in ref. [11] is likely to be more appropriate for the 
description of nonrelativistic particles in a rotating frame; cf. ref. [35]. 

Secondly, we have predicted the new GRE. This effect consists in the appearance of 
the electric current in the rotating matter composed of massless particles involved in the 
parity violating electroweak interaction. This electric current flows along the rotation axis. 
The new GRE is analogous to CME, known in QED, which consists in the generation of the 
electric current of massless charged particles along the external magnetic field [2, 3]. It should 
be noted that, for the first time, the analogy between the motion in a rotating electroweak 
matter and in an external magnetic field was mentioned in ref. [24]. 

Note that the appearance of the electric tension in a rotating conductor owing to the 
noninertial effects was also discussed in ref. [36]. The electric tension, predicted in ref. [36], 
is induced mainly by the Coriolis force acting on charged particles in a rotating conductor. 
If us is chosen along the z-axis, this tension is found in ref. [36] to be along the azimuthal 
direction. In our case, the electric current is owing to both the matter rotation and the 
presence of the parity violating electroweak interaction. We predict that the induced electric 
current flows along the rotation axis. 

We have used the solution of a Dirac equation in the rotating frame for the calculation of 
the induced electric current. It means that this current flows inside the rotating matter since 
the quantum states of charged particles are measured by a corotating observer. For example, 
if one used the wave functions obtained in ref. [37], although they look similar to those 
in eqs. (3.19) and (3.20), we would get the electric current with respect the a nonrotating 
observer, which cannot be applied for the generation of the internal TMF. 

We have used the calculated electric current to generate TMF and the magnetic helicity 
inside a rotating compact star, like NS or QS. The strength of TMF generated turned out 
to be moderate, B tor > 10 8 G, for both NS and QS. However, such TMF is comparable with 
a poloidal field in weakly magnetized old millisecond pulsars [20]. Note that the obtained 
strength of TMF is in agreement with a criterion for the magnetic field stability derived in 
ref. [18]. It should be noted that our model for the generation of TMF does not require the 
existence of a significant chiral imbalance between left and right charged particles. Such an 
imbalance is essential if CME is used to create TMF; cf. ref. [19]. 

Finally, in section 6, we analyzed the possibility to apply GRE to explain linear velocities 
of pulsars by the asymmetric neutrino emission from a rotating NS. We have estimated the 
total neutrino flux along the rotation axis as well as the recoil velocity of NS. It turned out 
that a pulsar kick caused by GRE is outside the observationally tested region. 

At the end of this section we should make a comment on the influence of nonzero 
masses of charged particles on the generation of an electric current in a rotating star. It was 
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mentioned in ref. [19] that the value of the current, induced by CME, slightly diminishes if a 
small, compared to the energy, but nonzero electron mass is accounted for. A more detailed 
analysis of the influence on the current, induced in frames of CME, from nonzero masses of 
flavored fermions, in case of strong interactions, was made in ref. [38]. It was shown in ref. [38] 
that CME disappears only in the great masses limit. The effect of the nonzero electron mass 
on the generation of strong magnetic fields in magnetars was also studied in ref. [39]. 

We should mention that, in case of GRE, a nonzero current in eq. (4.8) exists even at 
zero chiral imbalance: ps = 0. Therefore, even if an initial gts(0) 7 ^ 0 is washed out owing 
to spin-flip processes taking place at a nonzero mass, as predicted in refs. [4, 27, 39], the 
current will be nonzero due to V 5 = (El — Er)/ 2 / 0. For example, in section 5, we assumed 
that p§ = 0 to simplify the estimates. In the model for the generation of strong magnetic 
fields in magnetars, elaborated in refs. [4, 27], ~ —E 5 appears inside NS in the course of 

the magnetic fields evolution. In this case, our estimates for B t or obtained in section 5 will 
slightly change since gR 7 ^ /*l- 

GRE at V 5 / 0 and ^5 = 0 is likely to exist in all particular cases we considered in the 
present work since we discussed ultrarelativistic particles having rrif/(Ef) -C 1, where rrif 
is the fermion mass and ( Ef) is the typical fermion energy. Nevertheless, a more detailed 
quantum field theory analysis of this fact, like in ref. [38], is required. 
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A Details of the electric current calculation 


In this appendix we derive the electric current in eq. (4.8) and discuss the approximations 
made. 

Let us first consider left fermions. Using eq. (4.1), one gets that the contribution of 
particles and antiparticles to the electric current along the rotation axis is 



OO 


9/ 

N,s=0 


d p z 


^n tJ '(x)*Pfp 1 f(E f ) - -0/7 M (x)^/P/(-E/) , 


(A.l) 


where £/• j = £± Ll ± wsgn(VL )(2 ~ 1 / 2 ) are the energies of particles and antiparticles, which 
also include E c { = — (J • cj), and 0y j are the wave functions of particles and antipartic.les. 

Assuming that £ |E c f|, we can expand the distribution function in eq. (A.l) in a 
series 


Pf(Ef) = p L f (E 0 ) - 


wsgn(E L )Z' exp(£' 0 ) 


+ 


(A.2) 


T (exp(Eo) + l ) 2 

where Eq = £ ± El = Fwsgn(EL )/2 and V = l for the upper sign in Eq or l 1 = l — 1 for the lower 
sign in Eq. One can write down the decomposition of p^(Ej) analogous to that in eq. (A.2). 

Using the explicit form of the wave functions in eqs. (3.19) and (3.20) and the known 
sum involving the Laguerre function [40], 


J2(s-N)I 2 NjS (p) = p, 

s =0 


(A.3) 
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we get that accounting for E c f is equivalent to the replacement ^ —>• // L = + (u>r) 2 Vf. 
Here we neglected terms ~ w 3 . In fact, the terms ~ (cur) 2 are the noninertial contributions 
to the electric current. 

Thus the ground state N = 0 contribution to in eq. (A.l) takes the form, 


r3 _ „ wV L 

J L - -Qf 


IT 


J d p Pf(p + V L ) -p){p- VL) 


(A.4) 


where we should use // L instead of m,. Introducing the new variables x = p/T and y = 
(pL ~ Vl )/T in eq. (A.4) and using the identity 


J 0 (e x -y + 1 e x +y + l) dx y ’ 

and eq. (A.4), we get that J 3 = —{qf/ 7t)u;Vl(Ml — Vl). The contribution of right fermions 
can be obtained analogously. It is interesting to mention that the electric current obtained 
does not depend on the plasma temperature. 

In section 5 we use the following values of the parameters: pl,r ~ 10 2 MeV, Vl,r ~ 
10 eV, u) ~ 10 3 s _1 , and r < R ~ 10 km. Therefore, one gets that /tl,r Hl,r (wr) 2 I/L,R- 
Taking into account these estimates we arrive to eq. (4.8). 


References 

[1] L. Ryder, Introduction to General Relativity , Cambridge University Press, Cambridge U.K. 
(2009), pp. 18-46. 

[2] V.A. Miransky and I.A. Shovkovy, Quantum field theory in a magnetic field: From quantum 
chromodynamics to qraphene and Dirac semimetals, Phus. Rep. 576 (2015) 1, 

[arXiv: 1503.00732], 

[3] A. Vilenkin, Equilibrium parity-violating current in a magnetic field, Phys. Rev. D 22 (1980) 
3080. 

[4] M. Dvornikov and V.B. Semikoz, Magnetic field instability in a neutron star driven by the 
electroweak electron-nucleon interaction versus the chiral magnetic effect, Phys. Rev. D 91 
(2015) 061301 [arXiv: 1410.6676]. 

[5] A. Brandenburg and K. Subramanian, Astrophysical magnetic fields and nonlinear dynamo 
theory, Phys. Rept. 417 (2005) 1 [astro-ph/0405052]. 

[ 6 ] E. Flowers and M.A. Ruderman, Evolution of pulsar magnetic fields, Astrophys. J. 215 (1977) 
302. 

[7] R.J. Tayler, The adiabatic stability of stars containing magnetic fields-1. Toroidal fields, Mon. 
Not. R. Astron. Soc. 161 (1973) 365. 

[ 8 ] G. Hobbs, D.R. Lorimer, A.G. Lyne and M. Kramer, A statistical study of 233 pulsar proper 
motions, Mon. Not. R. Astron. Soc. 360 (2005) 974 [astro-ph/0504584]. 

[9] S. Johnston, G. Hobbs, S. Vigeland, M. Kramer, J.M. Weisberg and A.G. Lyne, Evidence for 
alignment of the rotation and velocity vectors in pulsars, Mon. Not. Roy. Astron. Soc. 364 
(2005) 1397 [astro-ph/0510260]. 

[10] R.N. Mohapatra and P.B. Pal, Massive Neutrinos in Physics and Astrophysics, World 
Scientific, Singapore (2004), 3rd. edn., pp. 5-8. 

[11] M. Dvornikov, Neutrino interaction with matter in a noninertial frame, JHEP 10 (2014) 053 
[arXiv: 1408.2735], 


- 15 - 







[12] L.D. Landau and E.M. Lifshitz, The Classical Theory of Fields, Butterworth Heinemann, 
Amsterdam (1994), 4th edn., pgs. 329-330. 

[13] D. Phiz, M. Roy and J. Wudka, Neutrino oscillations in strong gravitational fields, Phys. Rev. 
D 54 (1996) 1587 [hep-ph/9604403], 

[14] A.A. Grib, S.G. Mamaev and V.M. Mostepanenko, Quantum Effects in Intense External 
Fields: Methods and Results not Related to the Perturbation Theory, Atomizdat, Moscow 
U.S.S.R. (1980), pp. 13-15. 

[15] K. Bakke, Noninertial effects on the Dirac oscillator in a topological defect spacetime, Eur. 
Phys. J. Plus 127 (2012) 82 [arXiv:1209.0369], 

[16] P. Schluter, K.-H. Wietschorke and W. Greiner, The Dirac equation on orthogonal coordinate 
systems: I. The local representation, J. Phys. A: Math. Gen. 16 (1983) 1999. 

[17] C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill, New York U.S.A. (1980), 
pp. 691-696. 

[18] S.K. Lander and D.I. Jones, Magnetic fields in axisymmetric neutron stars, Mon. Not. R. 
Astron. Soc. 395 (2009) 2162 [arXiv:0903.0827]. 

[19] J. Charbonneau and A. Zhitnitsky, Topological currents in neutron stars: kicks, precession, 
toroidal fields, and magnetic helicity, JCAP 08 (2010) 010 [arXiv:0903.4450]. 

[20] E.S. Phinney and S.R. Kulkarni, Binary and millisecond pulsars, Annu. Rev. Astron. 
Astrophys. 32 (1994) 591. 

[211 F. Weber, Stranqe quark matter and compact stars, Proq. Part. Nucl. Phys. 54 (2005) 193 
[astro-ph/0407155]. 

[22] D.P. Menezes, C. Providencia and D.B. Melrose, Quark stars within relativistic models, J. 
Phys. G: Nucl. Part. Phys. 32 (2006) 1081 [astro-ph/0507529]. 

[23] C.Y. Ng, K.S. Cheng and M.C. Chu, Cooling properties of Cloudy Bag strange stars, 
Astropart. Phys. 19 (2003) 171 [astro-ph/0209016]. 

[24] A.V. Grigoriev, A.M. Savochkin and A.I. Studenikin, Quantum states of the neutrino in a 
nonuniformly moving medium, Russ. Phys. J. 50 (2007) 845. 

[25] M. Dvornikov, Impossibility of the strong magnetic fields generation in an electron-positron 
plasma, Phys. Rev. D 90 (2014) 041702 [arXiv:1405.3059], 

[26] M.A. Berger, Introduction to magnetic helicity, Plasma Phys. Control. Fusion 41 (1999) B167. 

[27] M. Dvornikov and V.B. Semikoz, Generation of the magnetic helicity in a neutron star driven 
by the electroweak electron-nucleon interaction, to be published in JCAP (2015) 

[arXiv: 1503.04162]. 

[28] E.R. Harrison and E. Tademaru, Acceleration of pulsars by asymmetric radiation, Astrophys. 
J. 201 (1975) 447. 

[29] A. Burrows, E. Livne, L. Dessart, C.D. Ott and J. Murphy, Features of the acoustic 
mechanism of core-collapse supernova explosions, Astrophys. J. 655 (2007) 416 
[astro-ph/0610175]. 

[30] A. Kusenko and G. Segre, Velocities of pulsars and neutrino oscillations, Phys. Rev. Lett. 77 
(1996) 4872 [hep-ph/9606428], 

[31] T. Totani, K. Sato, H.E. Dalhed and J.R. Wilson, Future detection of supernova neutrino 
burst and explosion mechanism, Astrophys. J. 496 (1998) 216 [astro-ph/9710203]. 

[32] M.T. Keil, G.G. Raffelt and H.-T. Janka, Monte Carlo study of supernova neutrino spectra 
formation, Astrophys. J. 590 (2003) 971 [astro-ph/0208035]. 


- 16 


[33] D.G. Yakovlev, W.C.G. Ho, P.S. Sternin, C.O. Heinke and A.Y. Potekhin, Cooling rates of 
neutron stars and the young neutron star in the Cassiopeia A supernova remnant , Mon. Not. 
Roy. Astron. Soc. 411 (2011) 1977 [arXive:1010.1154]. 

[34] F.W. Hehl and W.-T. Ni, Inertial effects of a Dirac particle , Phys. Rev. D 42 (1990) 2045. 

[35] K. Bakke and C. Furtado, Bound states for neutral particles in a rotating frame in the cosmic 
string spacetime , Phys. Rev. D 82 (2010) 2045. 

[36] B.J. Ahmedov, General relativistic galvano-gravitomagnetic effect in current carrying 
conductors , Phys. Lett. A 256 (1999) 9 [gr-qc/0701045]. 

[37] A.I. Studenikin and I.V. Tokarev, Millicharged neutrino with anomalous magnetic moment in 
rotating magnetized matter , Nucl. Phys. B 884 (2014) 396 [arXiv:1209.3245]. 

[381 C. Hoyos, T. Nishioka, and A. O’Bannon, A chiral maqnetic effect from AdS/CFT with flavor , 
JEEP 10 (2011) 084 [arXiv: 1106.4030], 

[39] D. Grabowska, D. Kaplan and S. Reddy, The role of the electron mass in damping chiral 
maqnetic instability in supernova and neutron stars, Phys. Rev. D 91 (2015) 085035 
[arXiv: 1409.3602], 

[40] A.A. Sokolov and I.M. Ternov, Relativistic Electron , Nauka, Moscow U.S.S.R. (1974), p. 354. 


- 17 - 


